A metamagnetic critical point in a three dimensional frustrated antiferromagnet 
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The competition between different forms of order is central to the problem of strong correla- 
tion. This is particularly true of frustrated systems, which frequently exist at or near to a zero- 
temperature critical point. Here we show that a state with a half-magnetization plateau but no long 
range order can arise when a three dimensional frustrated antiferromagnet is tuned to a critical point 
bordering a metamagnetic state. We use classical Monte Carlo simulation and low-temperature ex- 
pansion techniques to accurately characterize this "spin pseudogap" state, and show how its prop- 
erties relate to those of the critical point. Our results provide an example of three dimensional spin 
model which can be used to study the relationship between gap and "pseudogap" , — i.e. long range 
and preformed local order — near a metamagnetic critical point. 

PACS numbers: 75.10.-b, 75.10.Hk 75.80.+q 



In recent years the concept of "quantum criticality" 
has become central to efforts to understand correlated 
electron systems. This rapidly growing body of work 
rests on the simple idea that where there is a zero- 
temperature phase transition between two different or- 
dered phases, the finite temperature properties of the 
paramagnetic phase are controlled by the critical point 
which separates them. In order to test this hypothesis 
experimentally, it is necessary to identify systems with 
a suitable control parameter which can be used to tune 
through the critical point. Mechanical pressure, chem- 
ical pressure, doping and magnetic field have all been 
used successfully to this end. Two widely discussed 
examples are the underdoped cuprate superconductors, 
whose "pseudogap phase" has been suggested to origi- 
nate in a quantum critical point as a function of dop- 
ing, and Sr3Ru2 07 where non-Fermi liquid behavior ap- 
pears to be associated with a metamagnetic transition 
with strongly suppressed T c [lj. 

Theoretical attempts to understand these phenomena 
have largely concentrated on renormalization group anal- 
ysis of phenomenological field theories Q{. However it is 
also interesting to ask what happens in microscopic mod- 
els, especially where these are accessible to a variety of 
different approaches. Can we construct concrete exam- 
ples of systems with a zero-temperature critical point ? 
What are the nature of the correlations at this point, and 
in the paramagnetic phase connected to it ? What does 
a "pseudogap" look like at a microscopic level ? 

In this Letter we show how a state with a gap to spin 
excitations — a "spin pseudogap" — but no long-range 
magnetic order, can arise near a metamagnetic critical 
point in a microscopic model. Our main results are sum- 
marized in the phase diagram Fig.^ The model which we 
use to explore these ideas is the antiferromagnetic (AF) 
Heisenberg model on the highly frustrated pyrochlore lat- 
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FIG. 1: (Color online) Phase diagram for the classical Heisen- 
berg antiferromagnet on a pyrochlore lattice Q in applied 
magnetic field h = 4, with additional biquadratic interac- 
tions b — 0.6. The transition temperature T c associated with 
the gapped, ordered, half-magnetization plateau state van- 
ishes as the strength of ferromagnetic third-neighbour inter- 
actions J3 — » 0. A "spin-pseudogap" phase exhibiting a half- 
magnetization plateau but no long-range magnetic order ex- 
ists above this critical point below a crossover temperature 
T*. Solid lines are guides for Monte Carlo data, and the 
dashed line shows T c predicted by low-temperature expan- 



tice, in applied magnetic field h: 

(i,3) 

■ X Si-Si-h^S?. (1) 

In order to make the problem accessible to large scale 
simulation, we consider the classical limit of the problem 
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FIG. 2: (Color online) Half-magnetization plateau states 
(uuud states) on a pyrochlore lattice with exactly three up 
(black) and one down (white) spins per tetrahedron, (a) uuud 
state with long-range four-sublattice order, as considered in 
Q. (b) A schematic picture of uuud state with no long range 
order of the type found at the critical point. 




FIG. 3: (Color online) Dependence of (a) magnetization m 
and (b) uniform magnetic susceptibility \ on magnetic field 
h for b — 0.6 and J3 = 0, showing the existence of the magne- 
tization plateau in the absence of long-range magnetic order. 
Symbols (lines) denote the data for L = 16 (L = 8). The 
dashed lines show the T = results for long-range four- 
sublattice order. 



S = |S| — > 00. The dominant effect of quantum fluctua- 
tions in frustrated systems is known to be a tendency 
towards collinearity, and we characterize this through 
an effective biquadratic interaction b |j. An additional 
third-neighbour interaction J3 is used to tune the sys- 
tem away from a critical point at J3 = 0. The model JJJl 
was recently shown to offer an explanation of the broad 
half-magnetization plateaux observed in the spinel ox- 
ides CdCr 2 04 and HgCr 2 04 0,0- In these materials, 
strong effective biquadratic interactions b ~ 0.1-0.2 arise 
from the coupling of spins to the lattice. 

For simplicity, we focus below on the limiting case of 
ferromagnetic (FM) J 3 — >• 0~. For FM J 3 , the model (JTJ 
exhibits q = four-sublattice long-range magnetic order 
at low temperatures |4| . The most dramatic feature of its 
zero-temperature phase diagram is a first order transi- 
tion in applied field into a metamagnetic "plateau" state 
with m = 1/2 — exactly half the saturation magnetiza- 
tion — protected by a gap to the lowest lying spin excita- 
tion. This state is of the type illustrated in Fig. Eta). We 
have performed extensive Monte Carlo (MC) simulations 
of Eq. 0} at finite temperatures [(|, complemented by 
classical low-temperature expansions about both ordered 
and disordered states. These confirm that the phase di- 
agram given in Q remains valid up to a transition tem- 
perature T c whose energy scale is set by | J$\. However 
for I J3I — > 0, the transition temperature for each of the 
different ordered phases vanishes and the system exists 
at a critical point. 

So much for long range order — what about metamag- 
netism ? In Fig. we illustrate the magnetization process 
of Eq. (JTJ for J3 = 0. Data are taken from MC simula- 
tions of Eq. JTJ using a Metropolis algorithm with local 
spin update \J\. For convenience we assume a (large) 
value of b = 0.6, and work in units such that J = S = 1. 
It is clear from Fig. Ota) that the magnetization plateau 
is alive and well in fact from the magnetization pro- 
cess alone it is essentially impossible to distinguish these 



results from those for the four-sublattice ordered phase 
for FM J3 jfj. The strong suppression of the magnetic 
susceptibility at low temperatures shown in Fig.^b) sug- 
gests the existence of a well defined gap, of similar mag- 
nitude to that in the nearby ordered phase. How should 
we reconcile these results with the vanishing transition 
temperature for the ordered plateau state at J3 = ? 

In Fig. 2] we present MC results for the reduced spin- 
spin correlation function 



Q( rij ) = (Si ■ S, 
and the measure of collinearity 



((Si • Sj) 



(2) 



(3) 



As shown in Fig. Ufa) 7 spin correlations exhibit a liquid- 
like structure. After an initial AF oscillation on the scale 
of a near-neighbour they decay very rapidly to zero. 
Within the range accessible to MC simulation, the en- 
velope for this decay appears to cross over smoothly be- 
tween an exponential decay at high temperatures and a 
1/r 3 power-law behavior at low temperatures, as shown 
in Figs. 0tb) and (c). The collinearity P(ry), on the 
other hand, does exhibit long range order — see Figs.^Jd) 
and (e). 

We can understand these results as follows — the bi- 
quadratic interaction 6, which we have introduced to 
characterize quantum effects, favors collinear states. The 
pyrochlore lattice is corner-sharing network of tetrahe- 
dra. For h ~ 4, the b term acts to select uuud states 
with exactly three up and one down spins in each tetra- 
hedron and, consequently, m = 1/2. The q = ordered 
phase [Fig. Eta)] is an example of such states. But there 
are many more — in fact there is a one-to-one map- 
ping between uuud states on a pyrochlore lattice and the 
- 1.3 N / 2 dimcr coverings of the dual diamond lattice 81. 
Exactly at the critical point, for T = and J3 = 0, 
the system exists in an equally weighted superposition of 
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FIG. 4: (Color online) Absence of long-range magnetic or- 
der in the spin-pseudogap state for 6=0.6, J3 = and h—4. 
(a) — (c) The reduced spin correlation function Q(r), defined 
by Eq. 0, decays exponentially at high temperatures, cross- 
ing over to a power law behavior Q(r) ~ 1/r 3 at low tem- 
peratures, (d), (e) However the measure of collinearity P(r), 
defined by Eq. 10 does exhibit long range order at low tem- 
peratures. All the results are calculated along [flO] chains 
for L = 16. Distances are measured in units of the nearest- 
neighbor bond length. 



all uuud states, exhibiting a long-range spin collinear- 
ity but no magnetic order. A schematic picture of this 
spin-liquid plateau state is shown in Fig.Efb). 

We can use this insight into the nature of the ground 
state manifold at the critical point to construct a low- 
temperature theory for the paramagnetic phase con- 
nected with it. Expanding in small fluctuations about 
the dimer manifold, we find a free energy per spin 

T = J - 3bJ — h — TlnT + T(ln |M|)/2 

-(Tlnl.3)/2 + C(r 2 ), (4) 

where M is the Hamiltonian matrix for harmonic oscil- 
lations about a given uuud state, and the average (. . .) 
is taken in the manifold of all uuud states. Fluctua- 
tion effects, and in particular the magnetic susceptibility 
x(h,T), depend on the model parameters only through 
the entropy term T(ln |M|)/2. 

In the limit 6 — > 0, all uuud states have exactly the 
same excitation spectrum, with four distinct "bands" , 
two of them non-dispersing zero modes. Finite values 
of b lift the degeneracy of the spectra of different uuud 
states, and open a gap Ai to the first spin excitation. We 
have performed a MC search within the manifold of uuud 
configurations for finite values of b, evaluating the exci- 
tation spectrum and In |M| for each state numerically on 
a finite lattice of N = 1024 spins. Typical results for the 
cumulative density of states (DOS) are shown in Fig. 
We find that the spectrum is always bounded above and 
below by the same gaps A2 and Ai, which are set by 
the same nodeless excitations. In the four-sublattice or- 
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FIG. 5: (Color online) (a) DOS of four-sublattice uuud state 
in thermodynamic limit, showing finite gap and flat band at 
finite energy, (b) Cumulative DOS of four-sublattice uuud 
state in 1024 spin cluster, (c) Cumulative DOS of typical 
disordered uuud state in 1024 spin cluster. In all cases J3 = 0, 
h = 4, and b = 0.6. 



dered uuud state these reduce to the highest and lowest 
energy excitations with q = 0, respectively. One of the 
zero modes survives at finite b as a non-dispersing ex- 
citation at finite energy. Variation in entropy between 
different uuud states is negligible compared with the en- 
tropy of the uuud manifold, and so fluctuations alone 
cannot drive the system to order. The insensitivity of 
the excitation spectrum to the type of uuud state con- 
sidered explains why the thermodynamic properties of 
the disordered phase for J 3 = — in particular its half- 
magnetization plateau — are so similar to those of the 
ordered phases for finite | J3 1 . 

We can go further and calculate the asymptotic decay 
of spin correlations in the region of the critical point. 
We do this by constructing an effective electrodynamics 
for the uuud "dimer" manifold in close analogy with a 
recent treatment of the uudd "ice" manifold 9] . In this 
approach, the entropy term T In |M|/2 is written in terms 
of a course-grained polarization field P which measures 
deviations from uuud order. We find that the spin corre- 
lations in the uuud manifold decay as 1/r 3 , just like those 
of the uudd manifold. This result is in perfect agreement 
with the results of our MC simulation at low tempera- 
tures — see Fig. St c). 

Thus the situation in the immediate vicinity of the 
zero-temperature critical point for J3 = is fairly clear. 
A power law behavior of spin correlations coexists with a 
finite gap to spin excitations and a well defined plateau in 
the magnetization process. Since this gap is not associ- 
ated with long range order — at least in any conventional 
sense — we adopt the terminology "spin pseudogap" to 
describe this state. 

But how does this "spin pseudogap" state emerge from 
the paramagnet at high temperatures ? And how does it 
relate to the state for FM J3 with long-range magnetic 
order ? To answer these questions, we introduce mea- 
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FIG. 6: (Color online) Crossover into the spin-pseudogap 
phase as function of T for J3 = 0, h — 4 and b = 0.6 : Tem- 
perature dependence of (a) the specific heat (the gray line 
shows the contribution of the biquadratic term alone), (b) 
the local and global order parameters Aj 2 defined by Eq. 10 , 
and (c) the susceptibilities x\ associated with each of these 
order parameters. 



sures of both local and global order, based on the relevant 
irreducible representation of the tetrahedral symmetry 
group T d : 
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where, following |4|, Aj 2 = (A-r 2 ,i, A-r 2 ,2, At 2 ,3), and 
At s ,i = (S2 -S 3 — Si -Si)/V2, etc., are defined in terms of 
the spins on an individual tetrahedron. We also consider 
the related susceptibilities, \\{T) = [(A 2 ) - (A) 2 ] /T. 

In Fig. we show MC results for the specific heat, lo- 
cal and global "order" parameters Aj 2 and their related 
susceptibilities. A broad peak in the specific heat is ob- 
served at T ss b, with most of the entropy coming from 
the selection of collinear states. At the same time, the 
measure of local order Aj 2 cal increases rapidly, and the 
related susceptibility Xxt shows a broad peak. How- 
ever none of this signals the onset of long-range four- 
sublattice uuud order of the type found for FM J 3 < 
- the relevant order parameter Ay° a tends to zero in 
the thermodynamic limit, and the related susceptibility 
^global a pp earg 4- diverge only as T — » 0. 

The peak in x^ T cal or, equivalently, the peak in the 
specific heat c defines a crossover temperature T* for en- 
tering the spin pseudogap state, with its associated mag- 
netization plateau. The crossover also corresponds to an 
emergence of the power law behavior in spin correlations 
and the rapid development of collinearity in Fig. 0] T* 
is set by b — in fact, we have confirmed that for b < J, 
it scales nearly linearly with b. Meanwhile, the global 
order parameter monitors a transition for long-range or- 
dered phase for finite J3. By tracking both T* and T c 
as a function of J3 , we can map out the phase boundary 
of the true, long-range ordered plateau phase, and the 
extent of the finite-temperature spin pseudogap "phase" 
associated with the critical point at J3 = 0. Results for 
b = 0.6 and h — 4 are shown in Fig. ^ I n the vicinity 



of the critical point, we can also estimate T c for this (I s * 
order) transition within the low-temperature expansion 
by comparing the change in internal energy with the en- 
tropy change going from the gapped to the pseudogapped 
states: The prediction is T c ~ 30 1 J3 j, the dashed line in 
Fig. 2] This phase diagram summarizes the relationship 
between gap and "pseudogap" , or order and disorder, 
near a metamagnetic critical point. 

For small \J$\ where T c < T*, the system exhibits a 
magnetization plateau in the absence of long range or- 
der. As the system is cooled from the spin-pseudogap 
state into the ordered phase, the preformed local order 
in A-r 2 becomes global order in $Z tetra Aj 2 . This provides 
an intriguing magnetic analogy for the way in which pre- 
formed superconducting order is believed to emerge from 
preformed local cooper pairs in underdoped cuprates. 

So far as Cr spinels are concerned, our results suggest 
that a material which exhibits a magnetization plateau 
while maintaining the full (cubic) symmetry of the py- 
rochlore lattice, need not be magnetically ordered [lfj. 

In the light of our results it also seems highly proba- 
ble that the spin-1/2 Heisenberg model on a pyrochlore 
lattice exhibits a half-magnetization plateau because of 
strong quantum fluctuations. Within a spin-wave ap- 
proximation, our calculations show that the spin gap pro- 
tecting the plateau state is again finite and set by a node- 
less wave function, and so is insensitive to whether the 
system is ordered or disordered. 

We conclude by noting that it is possible to construct 
a phase diagram similar to Fig. for Eq. (JTJ near critical 
point with J3 = and h = 0. In this case the compet- 
ing ground states belong to the uudd "ice" manifold, and 
the correlations at the critical point are of the spin-liquid 
form described in 1 1 1 1 ] . However there are two important 
distinctions from the case of h = 4. Firstly, neither or- 
dered nor disordered uudd phases possesses a gap to spin 
excitations. And secondly, for h = 0, the collinearity 
Eq. |J3J is a well defined nematic order parameter. 

While we have been unable to construct any local ob- 
ject which can serve as order parameter for the spin- 
pseudogap state for h = 4, we cannot rule out the possi- 
bility that T* signals the onset of some exotic non-local 
order [T3 | . This remains as a subject for future investi- 
gation. 
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